Orbital-Free Density Functional Theory: Linear Scaling Methods 
for Kinetic Potentials, and Applications to Solid Al and Si 



o 
o 

(N 



o 



I 

O 

o 



> 
o 
in 

rn 
cn 

o 

On 
O 



X 



Jeng-Da Chai*,^ Vincent L. Ligneres,^ Gregory Ho,^ Emily A. Carter,'^ and John D. Weeks^'^ 

^Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742 
^Department of Chemistry, Princeton University, Princeton, New Jersey 08544 
"^Department of Mechanical and Aerospace Engineering and Program in Applied and Computational Mathematics, 

Princeton University, Princeton, New Jersey 08544 
^'Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742 

(Dated: March 20, 2009) 

In orbital- free density functional theory the kinetic potential (KP), the functional derivative of 
the kinetic energy density functional, appears in the Euler equation for the electron density and may 
be more amenable to simple approximations. We study properties of two solid-state systems, Al and 
Si, using two nonlocal KPs that gave good results for atoms. Very accurate results are found for Al, 
but results for Si are much less satisfactory, illustrating the general need for a better treatment of 
extended covalent systems. A different integration pathway in the KP formalism may prove useful 
in attacking this fundamental problem. 



I. INTRODUCTION 

Density-functional theory (DFT) is one of the most 
useful methods for studying ground state properties of 
many-electron systems. In principle, as shown by Ho- 
henberg and Kohn [l| , the exact ground state energy of a 
system of N electrons can be formally written as a func- 
tional E[p] of only the electron density p(r), a function 
of three variables, and the external field Vextiic) 0, 
Kohn and Sham (KS) [3| partitioned E[p] into the sum 
of three terms, the classical Hartree energy Eh[p] (readily 
expressed in terms of the density), the kinetic energy den- 
sity functional (KEDF) Tg [p] of a model system of non- 
interacting electrons at the same density, and the much 
smaller remainder, the so-called exchange-correlation en- 
ergy functional E^dp]- KS showed that the numerical 
value of Ts[p] could be determined exactly, not directly 
from the density itself, but by using a set of N one- 
electron wave functions (orbitals) to solve the N cou- 
pled KS equations that describe the noninteracting sys- 
tem. By using these results along with relatively simple 
approximations for the small -E2:c[p] term, the resulting 
KS-DFT has proven successful in many applications. 

However, because of the use of the KS orbitals, KS- 
DFT typically scales as 0{N^), which limits the ap- 
plicability of KS-DFT. By using localized orbitals some 
orbital-based linear-scaling methods have been developed 
that are especially successful for insulating systems with 
large band gaps @, Q. However, the prefactors of these 
types of approaches are usually so large that their com- 
putational costs become cheaper than the traditional KS- 
DFT (with 0{N^) scaling) only for very large systems. 

This problem could be circumvented if we could de- 
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velop an accurate "orbital-free" density functional theory 
(OF-DFT), where the kinetic energy is expressed in terms 
of the electron density alone [1, [1, H, [1] . Of course, given 
its large magnitude, this is a very difficult task, and in- 
deed simple local approximations like that used in the 
Thomas-Fermi (TF) model [H have proved very in- 
accurate. 

Through the recent efforts of many workers [ill, IHj 
Hi, m, M, [13, M, [H, M, [SH, [13, [H , we now have signif- 
icantly better approximations for the KEDF. There have 
been two main advances. The first is the use of nonlo- 
cal KEDFs that reproduce known exact results both for 
very slowly varying or very rapidly varying fields, and 
that give the exact linear response (LR) of the density of 
uniform model systems to small perturbations. The sec- 
ond is to focus on the more slowly-varying valence elec- 
tron density as described by weak local pseudopotentials 
[11 [11, [IE [13, [11, [11, [13, |3l|, [H m, m . For OF-DFT 
with LR-based KEDFs, the use of such pseudopotentials 
not only can reduce the computational cost relative to all- 
electron calculations, but also can improve its accuracy, 
since the valence system will be closer to the LR regime 
where Ts[p] is designed to be accurate 0, HI]. Very 
promising results using such OF-DFT methods have been 
obtained for a variety of nearly free-electron-like metals. 

However, a more accurate treatment of the KEDF is 
still needed in other applications with significant varia- 
tions in the density. The main problem is that the exact 
Ts[p] is highly nonlocal, and we have little idea of the 
functional form of the nonlocality for densities far from 
the linear response regime. 

In earlier work, Chai and Weeks (CW) [H, [sH] pro- 
posed a different way to attack this basic problem. In 
the usual approach, approximate forms for Ts[p] and 
Exc[p] are assumed, and the density p(r) is obtained from 
the variational principle (Euler equation) associated with 
minimizing E[p]: 

p^VTAr-dp]) + Veff{r;[p]). (1) 
The total energy of the inhomogeneous system is then 
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determined from the energy functional E[p\. All other 
physical quantities related to the ground-state density 
could also be computed. Here ^ is the chemical poten- 
tial (the Lagrange multiplier associated with the normal- 
ization condition J p(r)dr = N), and Ve//(r;[p]) is an 
effective one-body potential defined by 

K//(r;[p]) ^ Vff(r;[p]) + y.,(r;[p]) + K.t(r), (2) 

where 

VH{r;[p])^ I j^^dr' (3) 

is the Hartree potential, and T4c(r; [p]) = 6Ex:c[p]/Sp{r) is 
the exchange-correlation potential. Similarly we interpret 

VTAr;[p]) = ST,[p]/6p{r) (4) 

as the kinetic potential (KP) arising from the KEDF (bgI . 

113. 

In principle this process could be reversed, and the ex- 
act Ts [p] could be determined from the KP Vr^ (r; [p] ) 
by functional i nteg ration over density changes in all re- 
gions of space dJIH, m, [131 • Because of this integra- 
tion, Ts[p] is a functional that depends on the density 
everywhere even in the simplest case, exemplified by TF 
theory, where Vr, (r; [p]) is a simple function of the local 
density at r only. 

Of course, the true Vr^ (r; [p]) itself is much more com- 
plicated and must be a nonlocal functional of the density 
as well. However it seems reasonable to assume that it 
depends most strongly on the density in some local re- 
gion around r. Recent detailed arguments [H, |1§| sug- 
gest that the KP in this sense is more local than is the 
KEDF. Since most problems in devising accurate approx- 
imations for Ts[p\ have arisen from the nonlocality, this 
suggests that it could be worthwhile to try to develop 
approximations for the more local KP Vr^ (r; [p] ) . 

To that end, following ideas first introduced for the 
OF-KEDF's, CW proposed two different nonlocal KPs 
for atomic systems that reduce to known exact forms for 
both slowly varying and rapidly varying perturbations 
and also reproduce exact results for the linear response of 
the density of the homogeneous system to small perturba- 
tions [3^ . The latter requirement introduces nonlocality 
into the resulting KPs. CW also proposed two efhcicnt 
pathways to determine the numerical value of the kinetic 
energy by integration of the KPs, as discussed below in 

Sec. inn 

Using the CW KPs and the most accurate (Herring) 
pathway along with weak ab initio local pseudopoten- 
tials for the valence component of the electron density, 
CW obtained results for the total energies and ioniza- 
tion energies of atoms, and for the shell structure in the 
atomic radial density profiles that are in better agreement 
with calculations using the full Kohn-Sham theory than 
earlier KEDFs were able to produce. For example, the 
ionization energies of the first and second row atoms pre- 
dicted by the CW KPs have average errors of 1.8 eV and 



2.1 eV for the two different KPs, compared with the full 
KS-DFT, while the CAT KEDF [l^], a nonlocal KEDF, 
has an average error of 3.9 eV. Similarly, errors in the va- 
lence energies predicted by CAT KEDF are almost twice 
as large than those predicted by the CW KPs [34] . 

In this paper, we use the CW KPs to study two very 
different solid-state systems, a nearly-free-electron-like 
metal (solid Al) and a covalent material (solid Si). As 
discussed below, the original forms of these KPs use the 
local Fermi wave vector. This permits their use in atomic 
and molecular systems where the electron density van- 
ishes far from from the nuclei. However, this also leads 
to a quadratic scaling with system size [34] , which limits 
their usefulness for very large extended systems. Fol- 
lowing earlier work on the KEDF [l^, [3 , we introduce 
in Sec. [TT] an expansion method to reduce the compu- 
tational cost of evaluating these nonlocal KPs, and the 
truncated expansions essentially scale linearly with sys- 
tem size. The density pathway by which the kinetic en- 
ergy is obtained from the KP is described in Sec. IIIII 
Section IIVI compares the results of this method with the 
full KS-DFT and with other approaches using KEDFs 
and gives our conclusions. 

II. LINEAR SCALING METHODS 

CW introduced two different nonlocal KPs 34], one 
denoted LQ because it satisfies known exact results to 
second order at low wavevectors ("Low q"), and the other 
HQ with similar behavior at high wavevectors ("High 
q"). The nonlocality is generated by the requirement that 
these KPs also reproduce the exact linear response func- 
tion of a uniform electron gas at the the local Fermi wave 
vector (LFWV) kp{r), defined in terms of the electron 
density by 

kp{r) EE {3n'p{r)Y/\ (5) 

Both potentials can be compactly described using the 
following generalized form with different values for a pa- 
rameter a: 

K(r; [p],kF{r)) = VrFir; [p]) + Vw{r; [p]) 

+ ^CFpP{v) j f{\v~v'\-kF{v))p"{v')dv' 

Here (3 = (2/3 — a). The HQ model uses a = 2/3 and the 
LQ model has a — 1/2. Vtf{y\ [p]) and Vwiv; [p]) are the 
TF KP and the von Weizsacker (W) KP [4^ respectively, 
and Cf is a numerical coefficient associated with the TF 
KEDF. See Ref. fsll for a detailed discussion. 

The weight function /(|r — r'|; A;i^(r)) is defined by the 
following inverse Fourier-transform-like integral: 

f{\v-v'\:kF{v))^-^ J />/2fc;^(r))e-*-('-^')dk. 

(7) 

Here 

fiq) = FLiq)-3q'-l. (8) 
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is directly related to the inverse linear response function 
F]^{q) of the uniform electron gas with density po re- 
duced wave vector q = k/2kpo, with kpo = i'^'^^PoY^^ 
El: 



FUq) 



In 



1 



1 



(9) 



Clearly kpo is replaced by the LFWV in the the defi- 
nition of the weight function /(|r — r'|; ^^(r)) in Eq. ([7]). 
Given this definition it is straightforward to rewrite the 
last term in Eq. ([6]) in Fourier space as 
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Cpp^ir) / /(fc/2fcf(r))p"(k)e-^'''-dk. (10) 



Because of the ^^(r) term in Eq. pi}|) . the integration 
must be done on a grid over a range of r, and all these 
models scale quadratically. 

In order to achieve the desirable linear scaling for fast 
computation in large solid-state systems, we follow the 
work of Wang et al. [l^ and apply a Taylor series ex- 
pansion to the weight functional /(|r — r'|;fci?(r)) with 
respect to a reference density p*. For extended systems, 
where p(r) is not significantly different from the average 
density po, the natural choice of is po. Throughout 
this paper, we choose p* — po in all of our calculations. 

The weight function /(|r — r'|; fci?(r)) can be expanded 
with respect to p* as: 



/(|r-r'|;fcf(r)) = /(|r-r'|;fc» 

+/(i)(|r-r'|;fc|,)(p(r)-p,) 
+ i/W(|r-r'|;fcJ,)(p(r)-p,)2 + 



(11) 



where k*p = {iir'^ p^f^ , and /(") (|r - r'|; A:^) = 9"/(|r- 
r'l; fci?(r))/i9p"(r) is the nth order derivative of /(|r — 
r'|;fci?(r)) with respect to p(r), and is evaluated at p*. 
Their functional forms, up to second order, in reciprocal 
space are 



q* 
3p, 



and 



qlf"{q*) + 4q*f'{q*) 
(3p.)2 



(12) 



(13) 



where = k/{2kp), and f'{q*) and f"{q*) are the first 
and the second derivative of /(g) with respect to g, . Since 
the analytical form of f{q) is available in Eq. ([8]), all the 
terms needed in the Taylor series expansion f^^^\q) can 
be obtained analytically. This simplicity is one of the 
main advantages of the KP method over related methods 
that use the KEDF, where these terms have usually been 
be obtained numerically, with one recent exception [4^ . 

By carrying out the Taylor series expansion to nth or- 
der in Eq. (|lip . and inserting it in Eq. Q, one can then 
take out the (p(r) — p*)™ factor of the mth term from its 



integral (where m is a nonnegative integer, and m < n), 
and the remaining integral becomes: 



/('")(|r-r'|;fc» p"(r')dr' 



(14) 



The integral in Eq. p4)l can be easily evaluated by a sim- 
ple fast Fourier transform (FFT). Therefore, to compute 
the nonlocal term in Eq. ^ using this scheme, one needs 
to evaluate a total of n -f 2 FFT's, including the FFT of 
p". Since all of the terms can be computed by FFT's, this 
scheme essentially scales linearly 0{M h\M) with system 
size, where M is the number of grid points. 

The corresponding linear-scaling HQ {a = 2/3) and 
LQ (a = 1/2) KPs with p, = po can be written as: 

K^™(r; [p], /cj.o) = Vtf{v. [p]) + VW(r; [p]) 
+ ^CppP{v) I j /(|r - r'l; fc^^o) p"(r')dr' 

-f (p(r) - po) / f'\\v - r'l; fc^^o) p"(r')dr' 



+ i(p(r) - po)2 / /(2)(|r - r'l; fc^^o) p"(r')rfr' + 



(15) 

where kpQ = (37r^po)^''^ is the uniform Fermi wave vec- 
tor. For simple nearly-free-electron metals like Al, the 
expansion method works very well, and indeed often 
only the zeroth order term is needed. For extended sys- 
tems with large density variations over space, this ex- 
pansion method could experience convergence problems 
similar to those that have been seen by other workers 

[13, EH a 111. 

However, in those cases, both the use of the LFWV 
in Eq. [5] and the basic linear response treatment of the 
nonlocality are probably inadequate as well. A two-body 
Fermi wavevector kpir, r') like the one used in the WGC 
KEDF 112] could be introduced, but this functional form 
represents an additional approximation that does not sys- 
tematically improve the underlying linear response treat- 
ment of the nonlocality. As we will see, the treatment of 
extended systems with large density variations remains a 
major challenge for all OF-DFT methods. 



III. DENSITY PATHWAY 

As discussed in Ref. [11], an integration pathway is 
needed to determine the value of Ts[p\ from a given 
Vt, (r; [p]) [H, 0, m, E^. If the exact Vt, (r; [p]) is used 
and the integration is carried out exactly, then all path- 
ways would give the same exact result for Tg [p] . An ap- 
proximate KP can give the same (approximate) value 
for the kinetic energy independent of pathway only if 
it arises from functionally differentiating a KEDF (the 
usual OFDFT approach) or equivalently if it exactly sat- 
isfies the nonlocal consistency conditions given by Her- 
ring [23 in his Eq. (24). As he points out, it is not easy 
to satisfy these formal conditions for general functionals, 
and we do not try to do so here. 
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Instead, we generate nonlocality by imposing a more 
easily implemented and physically suggestive condition, 
requiring that linear response theory is exactly satisfied 
for small perturbations about a uniform system. Thus 
our results will be independent of pathway in the linear 
regime but the kinetic energy obtained from integration 
of our approximate KPs for systems with large density 
variations will in general depend on the particular inte- 
gration pathway used, as discussed in detail in Ref. [33 |. 
Experience with the classical analogue of the KP method 
applied to nonuniform hard sphere fluids [4^ has shown 
that errors from this path dependence can be small (e.g., 
less than one percent for the surface tension at a hard 
wall) when used with particular pathways that do not 
excessively weight regions poorly described by an approx- 
imate theory for the nonuniform density. 

Two efficient pathways to obtain the kinetic energy 
from an approximate KP for atomic and molecular sys- 
tems were discussed in Ref. i34 | . The simplest and most 
accurate pathway, due to Herring (22| , automatically sat- 
isfies the virial theorem. However it does not apply to ex- 
tended systems in its original form and we have not yet 
found the appropriate generalization. Here we use the 
alternate linear density pathway, where there is a linear 
scaling of the electron density 



Table I: Lattice parameters (A) for bulk Al. The results for 
the KS, WT, and WGC models are taken from Ref. fTl. 



px{r) = PQ + X[p(r) - po] 



(16) 



If the exact KP were used, both pathways would give the 
exact kinetic energy. 

By inserting the HQ and LQ KPs of Eq. (fT5| into the 
linear density pathway (see Eqs. (9)-(ll) in Ref. [11]), 
the kinetic energy can be computed as 

T„ = Ta=o + Jo dXjdr[p{v) ~ poM"'{r; [px],kFo). 

(17) 

Here Tx=o is the kinetic energy for uniform system, i.e., 
the Thomas- Fermi kinetic energy Ttf[po]- Since the 
VTF(r; [px]) and Vw(r; [p\]) in Eq. p?)) arise from the 
functional derivatives of the known Ttf[p] and IV [p] 
functionals respectively, these terms can be integrated 
exactly and lead to the TF and W KEDFs. 
Equation pT]) then becomes 

Ta[p] = TtfIp] + Tw[p] + dxj rfr[p(r) - po] 
x^CFpiir)y dr'fi\r-r'\;kFo)ptir') 
+ {p,{r) - po) Jdr'/f^dr - r'\; kFo)pt{r') 

+ i(p,(r) - Po)' I dr'f('\\r - r'|; fc^o)pX(r') + ■ • j 

(18) 

where T\=o — Ttf[po] term is absorbed in the Ttf\p\ 
term. We use Eq. to evaluate the kinetic energy of 
our linear-scaling HQ (a — 2/3) and LQ [a = 1/2) KPs 
in all the calculations here. 



Al 


KS 


LQ 


HQ 


WT 


WGC 


fee 


4.03 


4.04 


4.04 


4.04 


4.03 


bee 


3.23 


3.23 


3.23 


3.23 


3.22 


sc 


5.33 


5.31 


5.36 


5.33 


5.38 


dia 


5.84 


5.91 


5.90 


5.94 


5.92 



Table II: Energy per atom (eV) for bulk Al. The first row is 
the energy for the fee structure, while other rows are energy 
difference from the fee structure. The results for the KS, WT, 
and WGC models are taken from Ref. 1131. 



Al 


KS 


LQ 


HQ 


WT 


WGC 


fee 


-58.336 


-58.303 


-58.314 


-58.331 


-58.331 


bee 


0.068 


0.053 


0.057 


0.060 


0.066 


sc 


0.250 


0.253 


0.253 


0.227 


0.217 


dia 


0.599 


0.712 


0.751 


0.673 


0.584 



IV. RESULTS 

In solids, the external potential Vext{^) can be regarded 
as a linear combination of the special array of local atomic 
pseudopotentials centered at each ion position R/. Dif- 
ferent arrays of R/ lead to different phases, such as face- 
centered cubic (fee), diamond (dia), body-centered cu- 
bic (bcc), simple cubic (sc), and so on. We calculate 
the binding energies of Al and Si at these four different 
phases, and compare our results with other KEDFs and 
KS-DFT. We use periodic boundary conditions with a 
cubic supercell containing 4 atoms for fee, 8 atoms for 
dia, 2 atoms for bcc, and 8 atoms for sc. All calculations 
are spin-restricted and use the local density approxima- 
tion (LDA) [13, m, m, [131 for the exchange-correlation 
functional. 

For bulk Al, the empirical Goodwin-Needs-Heine 
(GNH) local pseudopotential [l^l is used, and a plane 
wave kinetic energy cutoff of 600 eV is used to converge 
the electron density. Here it is sufficient to use only the 
zeroth-order linear-scaling HQ and LQ KPs (sec Eq. ^TS]) ) 
in the calculations, as verified by comparison with results 
of a first order calculation. Our results are compared with 
the WT [H and WGC [iIkEDFs, and KS-DFT, which 
were previously computed |12l | , using the same local pseu- 
dopotential. (The atomic-based local pseudopotentials in 
Ref. [13| could also be used and gives very similar results.) 
As can be seen in Tables HI and iHl all the LR-based mod- 
els perform similarly, and agree well with KS-DFT. The 
phase ordering is correct, and aside from the high energy 
diamond phase, the lattice parameters are close to the 
KS results. 

To assess the performance of the HQ and LQ KPs in 
covalent systems like solid Si, we use the bulk local pseu- 



dopotential (BLPS) developed by Zhou et al. [2j| to- 
gether with a plane wave kinetic energy cutoff of 2000 
eV for converging the electron density. We used the first- 
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Table III: Lattice parameters (A) for bulk Si. The results for 
the KS and WGC models are taken from Ref. [2J, and those 
for the WGC2 model are taken from Ref. fTi|. 



Si 


KS 


LQ 


HQ 


WGC 


WGC2 


dia 


5.38 


5.27 


5.29 


5.77 


5.57 


bcc 


3.29 


3.06 


3.06 


3.29 


3.32 


sc 


4.99 


4.98 


4.98 


5.01 


5.06 


fee 


3.83 


3.83 


3.82 


3.80 


3.80 



Table IV: Energy per atom (eV) for bulk Si. The first row 
is the energy for the dia structure, while other rows are en- 
ergy difference from the dia structure. The results for the KS 
and WGC models are taken from Ref. 24], and those for the 
WGC2 model are taken from Ref. 14] . 



Si 


KS 


LQ 


HQ 


WGC 


WGC2 


dia 


-110.234 


-109.167 


-109.282 


-110.345 


-110.220 


bcc 


0.165 


-0.553 


-0.437 


0.537 


0.267 


sc 


0.303 


-0.586 


-0.531 


0.506 


0.308 


fee 


0.457 


-0.584 


-0.478 


0.571 


0.444 



order linear-scaling HQ and LQ KPs in the Si calcula- 
tions. Our results are compared first with the original 
WGC KEDF [13, and with the KS-DFT, which were 
previously computed [IJ] , using the same local pseudopo- 
tential. As can be seen in Tables Hill and [TVl all these LR- 
based models perform worse than for Al, when compared 
with KS-DFT. The phase ordering is incorrect, and the 
lattice parameters of the four phases only qualitatively 
match with the KS results, although the WGC KEDF at 
least obtains the correct diamond ground state. 

Because we use approximate KPs, our estimate for 
Ts[p] will depend on the integration pathway, and it is 
important to try to find a particular pathway that is rel- 
atively insensitive to the errors that exist in our KPs. 
For atomic systems, we were able to use the Herring 
pathway, which automatically satisfies the virial theorem, 
and good results were found [34]. As noted above, a gen- 
eralization of the Herring pathway for solids is not yet 
available, and in this paper we used the linear density 
pathway. While this pathway is numerically efficient, 
and gives exact results when the exact KP is used, we 
expect less accurate results when used with our approxi- 
mate KPs since the virial theorem is not satisfied. 

Indeed this is the case for atomic systems, where results 
from both pathways can be compared. In Tables |Vl we 
show the pathway dependence of the full LQ KP (with 
ab initio local pseudopotentials [s^ D in atomic systems. 
The kinetic energies of the three atoms (Al, Si and Ar), 
using the full LQ KP, are evaluated by the Herring and 
the density pathways, and are compared with the KS 
results. Clearly, the density pathway gives significantly 
worse kinetic energy than the Herring's pathway, and its 
errors increase when the densities of the systems are more 
rapidly varying. This suggests that the large deviations 
of the LQ and HQ results for solid Si (dia) in Tables IIVI 
may be partially due to the use of this pathway. This 



Table V: The kinetic energy T^fp] (eV) of the three atoms 
(Al , Si, and Ar) in ab initio local pseudopotential calculations 
[34| . using the KS method and the full LQ model (evaluated 
by both the Herring's and the density pathways). MAE, the 
mean absolute errors (relative to the KS method) of the full 
LQ model computed by the two pathways are given at the 
bottom of their respective columns. 





KS 


LQh 




Al 


21.164 


21.823 


17.563 


Si 


40.985 


41.264 


33.874 


Ar 


225.609 


221.035 


196.495 


MAE 




1.838 


13.276 



pathway dependence also seems likely to influence the 
phase orderings. Thus the generalization of the accurate 
Herring pathway to solids seems a promising focus for 
future research. 

Since Si is a covalent material, the density inside the 
covalent bond regions is quite different from that out- 
side. Even with a local pseudopotential this system may 
be outside the linear response regime, or the results at 
least may depend sensitively of the choice of the one- or 
two-body FWV. In particular, recent work 14] has shown 
that much improved results for the phase energies can be 
obtained for bulk Si using the WGC KEDF with a FWV 
mixing parameter and reference density optimized for the 
covalent phases of Si, rather than the original FWV mix- 
ing parameter optimized for Al p^ . Results with this 
new parameter set, denoted here WGC2, are also given in 
Tables IIIII and llVi where we see that the phase orderings 
and energy differences for Si phases using WGC2 are in 
very good agreement with KS-DFT. However, ultimately 
it is desirable to find a general KEDF or KP that does not 
depend on special properties of the system, particularly 
when considering more complicated systems with defects 
or surfaces. In our opinion, developing a truly univer- 
sal and yet accurate KEDF or KP for covalent materials 
remains an outstanding problem. 



V. CONCLUSION 

In summary, CW have previously demonstrated that 
the nonlocal HQ and LQ KPs work well for isolated atoms 
and ions [34| , and now we show they can be used with no 
additional parameterization in solid-state systems. Very 
good results are found for Al. However, problems of inac- 
curate phase ordering and bond lengths in Si are found, 
and a better treatment of nonlocality beyond the linear 
response regime may be needed for covalent systems with 
significant density variations. The simpler and more lo- 
cal form of the OF KP could be useful in future devel- 
opments. If it is possible to develop a new pathway to 
obtain the kinetic energy in extended systems that sat- 
isfies the virial theorem, analogous to the accurate Her- 
ring pathway for atomic systems [3^ , we believe even the 
present LR-based KP's would likely give more accurate 
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phase energies. Further work along these lines is called 
for. 
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